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Recently, there has been renewed interest in the search for low-inass magnetic 

monopoles. At the University of Oklahoma we are performing an experiment 
OO \ (Fermilab E882) using material from the old DO and CDF detectors to set 

' limits on the existence of Dirac monopoles of masses of the order of 500 GeV. 

To set such limits, estimates must be made of the production rate of such 
, J^ ■ monopoles at the Tevatron collider, and of the binding strength of any such 

^S| ■ produced monopoles to matter. Here we sketch the still primitive theory 

Y~l . of such interactions, and indicate why we believe a credible limit may still 

f^ \ be obtained. On the other hand, there have been proposals that the classic 

On ' Euler-Heisenberg Lagrangian together with duality could be employed to set 

limits on magnetic monopoles having masses less than 1 TeV, based on virtual. 






rather than real processes. The DO collaboration at Fermilab has used such 
a proposal to set mass limits based on the nonobservation of pairs of photons 



Q , each with high transverse momentum. We critique the underlying theory, 

'-^ \ by showing that the cross section violates unitarity at the quoted limits and 

is unstable with respect to radiative corrections. We therefore believe that 
^ ■ no significant limit can be obtained from the current experiments, based on 

H , virtual monopole processes. 
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I. INTRODUCTION 

The notion of magnetic charge has intrigued physicists since Dirac [^ showed that it was 
consistent with quantum mechanics provided a suitable quantization condition was satisfied: 
For a monopole of magnetic charge g in the presence of an electric charge e, that quantization 
condition is (in this paper we use rationalized units) 

where N is an integer. For a pair of dyons, that is, particles carrying both electric and 
magnetic charge, the quantization condition is replaced by 0,0 

An = y^^' (2) 

where {ei, gi) and (62,(72) are the charges of the two dyons.[] 

With the advent of "more unified" non-Abelian theories, classical composite monopole 
solutions were discovered 0]. The mass of these monopoles would be of the order of the 
relevant gauge-symmetry breaking scale, which for grand unified theories is of order 10^^ 
GeV or higher. But there are models where the electroweak symmetry breaking can give 
rise to monopoles of mass ~ 10 TeV |0. Even the latter are not yet accessible to accelerator 
experiments, so limits on heavy monopoles depend either on cosmological considerations 
0, or detection of cosmologically produced (relic) monopoles impinging upon the earth or 
moon 0. However, a priori, there is no reason that Dirac/Schwinger monopoles or dyons 
of arbitrary mass might not exist: In this respect, it is important to set limits below the 1 
TeV scale. 

Such an experiment is currently in progress at the University of Oklahoma p], where 
we expect to be able to set limits on direct monopole production at Fermilab up to several 
hundred GeV. This will be a substantial improvement over previous limits |^. But indirect 
searches have been proposed and carried out as well. De Riijula ||10[ proposed looking 
at the three-photon decay of the Z boson, where the process proceeds through a virtual 
monopole loop. If we use his formula |T^ for the branching ratio for the Z — > 87 process. 



compared to the current experimental upper limit [^ for the branching ratio of 10~^, we 
can rule out monopole masses lower than about 400 GeV, rather than the 600 GeV quoted 
in Ref. [jlO|- Similarly, Ginzburg and Panfil [jl2| and more recently Ginzburg and Schiller 



P^ considered the production of two photons with high transverse momenta by the collision 
of two photons produced either from e~^e~ or quark- (anti-) quark collisions. Again the final 
photons are produced through a virtual monopole loop. Based on this theoretical scheme. 



an experimental limit has appeared by the DO collaboration [|IJ], which sets the following 
bounds on the monopole mass M: 



^An additional factor of 2 appears on the right hand side of these conditions if a symmetrical 
solution is adopted — see Eq. (Il6|) below. 



^ r 610 GeV for ^ = 

— ><^ 870 GeV for 5 = 1/2 , (3) 

^ [ 1580 GeV for 5 = 1 

where S is the spin of the monopole. It is worth noting that a mass hmit of 120 GeV for 



a Dirac monopole has been set by Graf, Schafer, and Greiner [T^, based on the monopole 
contribution to the vacuum polarization correction to the muon anomalous magnetic mo- 
ment. (Actually, we believe that the correct limit, obtained from the well-known textbook 
formula |jl6| for the gf-factor correction due to a massive Dirac particle is 60 GeV.) 



The purpose of the present paper is to, first, describe the key theoretical elements nec- 
essary for the establishment of a direct limit for production of monopoles at Fermilab: An 
estimate of the production cross section for monopole-antimonopole pairs at the collider 
must be made, and then an estimate of the binding probability of such produced monopoles 
with matter must be given, so that we can predict how many monopoles would be bound 
to the detector elements that are run through our induction detector. Such estimates were 
given in our proposal to Fermilab; the complete analysis will be given in the experimental 
papers to follow. Here the emphasis will be on the elementary processes involved. 



A secondary purpose of this paper is to critique the theory of Refs. |T^, ]12[, [|T3l, and 



T5| , and thereby demonstrate that experimental limits, such as that of Ref. [|1^], based on 
virtual processes are unreliable. We will show that the theory is based on a naive application 
of electromagnetic duality; the resulting cross section cannot be valid because unitarity is 
violated for monopole masses as low as the quoted limits, and the process is subject to 



enormous, uncontrollable radiative corrections. It is not correct, in any sense, as Refs. |T^ 
and 10 state, that the effective expansion parameter is gu/M, where u is some external 
photon energy; rather, the factors of uj/M emerge kinematically from the requirements of 
gauge invariance at the one-loop level. If, in fact, a correct calculation introduced such 
additional factors of u/M, arising from the complicated coupling of magnetic charge to 
photons, we argue that no limit could be deduced for monopole masses from the current 
experiments. It may even be the case, based on preliminary field-theoretic calculations, that 
processes involving the production of real photons vanish. 

II. EIKONAL APPROXIMATION FOR ELECTRON- (OR QUARK-) 

MONOPOLE SCATTERING 

It is envisaged that if monopoles are sufficiently hght, they would be produced by a 
Drell-Yan type of process occurring in pp collisions at the Tevatron. The difficulty is to 
make a believable estimate of the elementary process gg — > 7* — > MM, where q stands 
for quark and M for magnetic monopole. It is not known how to calculate such a process 
using perturbation theory; indeed, perturbation theory is inapplicable to monopole processes 
because of the quantization condition (|l|). It is only because of that consistency condition 
that the Dirac string, for example, disappears from the result. 

Only formally has it been shown that the quantum field theory of electric and magnetic 
charges is independent of the string orientation, or, more generally, is gauge and rotationally 
invariant P,|l3- It has not yet proved possible to develop generally consistent schemes for 



calculating processes involving real or virtual magnetically charged particles. Partly this is 
because a sufficiently general field theoretic formulation has not yet been given; this defect 
will be addressed elsewhere |T^. However, the nonrelativistic scattering of magnetically 
charged particles is well understood^ [|^,|^. Thus it should not be surprising that an 
eikonal approximation gives a string-independent result for electron-monopole scattering 
provided the condition (^ is satisfied. More than two decades ago Schwinger proposed 
0] and Urrutia carried out such a calculation |]2T[]. Indeed, this subject has had arrested 



development. Since this is the only successful field-theoretic calculation yet presented, it 
may be useful to review it here. (A more detailed discussion will be presented in Ref. [0.) 
The interaction between electric (J^) and magnetic (* J^) currents is given by 

l^M = -e^^^r f{dx){dx'){dx")J^'{x)&'D+{x - x')r{x' - x"yr{x"). (4) 



Here D+ is the usual photon propagator, and the arbitrary "string" function f^{x — x') 
satisfies 

d^f''{x-x') = 6{x-x'). (5) 

It turns out to be convenient for this calculation to choose a symmetrical string, which 
satisfies 

Fix) = -n-x). (6) 

In the following we choose a string lying along the straight line n^, in which case the function 
may be written as a Fourier transform 

rv r m_ ,^ / 1 1 \ 

•^'^^ ^ 2U (27r)4'^ {n-k-ie^ n-k + tej- ^^ 

In the high-energy, low-momentum-transfer regime, the scattering amplitude between 
electron and monopole is obtained from Eq. (H) by inserting the classical currents. 



oo ^/^ 



r(x) = el dX — 6(x-—x], (8a) 

J-oo m \ m J 

where m and M are the masses of the electron and monopole, respectively. Let us choose a 
coordinate system such that the incident momenta of the two particles have spatial compo- 
nents along the z-axis: 

P2 = (p, 0,0,p), P2 = (p,0,0, -p), (9) 



^Contrary to the statement in the second reference in Ref. |jl^, the nonrelativistic calculation is 
exact, employs the quantization condition, and uses no "unjustified extra prescription." 



and the impact parameter lies in the xy plane: 

6=(0,b,0). (10) 

Apart from kinematical factors, the scattering amplitude is simply the transverse Fourier 
transform of the eikonal phase, 

/(q) = I £h e-'^"^ (e'^ - l) , (11) 

where x is simply W'^'^^'^ with the classical currents substituted, and q is the momentum 
transfer. 

First we calculate X'-, it is immediately seen to be, if n^ has no time component, 

eg f d^k^ z • (n x k^) ^^^^.^ ^1,1 



^ 2 7(277)2 kl-te ""' ^ [n-k^-te^ h-k^ + iej' ^^^^ 

where k^ is the component of the photon momentum perpendicular to the z axis. From this 
expression we see that the result is independent of the angle n makes with the z axis. We 
next use proper-time representations for the denominators in Eq. (p!2D, 

X /"OO 



k 



2 







t/se-*^i, (13a) 



1 1 1 

-zr—, - + ^^-i = - 

n ■ k^ — 2e n ■ k_|_ + «e t 



/■O 

dX — dX 

"'-OO 



g«An.kxg-|A|e_ ^^^g^^^ 



We then complete the square in the exponential and perform the Gaussian integration to 
obtain 

or 

^^'JLt,n-^i^A±\. (15) 

27r \z ■ (b X n)y 

Because e*^ must be continuous when n and b lie in the same direction, we must have the 
Schwinger quantization condition for an infinite string, 

eg = AnN, (16) 

where A^ is an integer. 

To carry out the integration in Eq. ([TTD , choose b to make an angle ip with q, and the 
projection of n in the xy plane to make an angle with q; then 

X = |f(^- 0-^/2). (17) 

To avoid the appearance of a Bessel function, we first integrate over b = |b|, and then over 



/(q) = / d^ Jj^Jj^-ibq{cos^-ie)^2iN{i,-4>-n/2) 

Jo Jo 

4 g-2JAf{0+7r/2) ^ ^^^2N~1 



i q^ Jc {z + l/z — ieY 

-^---^ (18) 



g2 



where C is a unit circle about the origin, and where again the quantization condition (jTBD 
has been used. Squaring this and putting in the kinematical factors we obtain Urrutia's 
result (gg 



da {egf 1 ^ ^ 



dt 4tt f^' 



q\ (19) 



which is exactly the same as the nonrelativistic, small angle result found, for example, 
in Ref. |jl9|. This calculation, however, points the way toward a proper relativistic treat- 



ment, and will be extended to the crossed process, quark- antiquark production of monopole- 
antimonopole pairs elsewhere. 

III. BINDING OF MONOPOLES TO MATTER 

Once the monopoles are produced in a collision at the Tevatron, they travel through the 
detector, losing energy in a well-known manner (see, e.g., Ref. [^), presumably ranging out. 



and eventually binding to matter in the detector (Be, Al, Pb, for example). The purpose of 
this section is to review the theory of the binding of magnetic charges to matter. 

We consider the binding of a monopole of magnetic charge g io a, nucleus of charge Ze, 
mass M. = Arrip, and magnetic moment 

^^ = — 7S, (20) 

nip 

S being the spin of the nucleus. (We will assume here that the monopole mass ^ Ai, which 
restriction could be easily removed.) Other notations for the magnetic moment are 

7 = l + /. = f. (21) 

The charge quantization condition is given by Eq. (|lD. Because the nuclear charge is Ze, 
the relevant angular momentum quantum number is [recall A^ is the magnetic charge quan- 
tization number in Eq. (|1])] 

I - ^. (22) 

We do not address the issue of dyons [^], which for the correct sign of the electric charge 
will always bind electrically to nuclei. 



A. Nonrelativistic binding for S = 1/2 

In this subsection we follow the early work of Malkus 1^^ and the more recent paper of 



Bracci and Fiorentini p5|]. (There are also the results given in Ref. p6|, but this reference 
seems to contain errors.) 

The neutron {Z = 0) is a special case. Binding will occur in the lowest angular momen- 
tum state, J = 1/2, if 



l7l > 



2N 



(23) 



Since 7„ = —1.91, this condition is satisfied for all A^. 

In general, it is convenient to define a reduced gyromagnetic ratio, 

. A ^ 

7 = — 7, K = 7-l. 



(24) 



This expresses the magnetic moment in terms of the mass and charge of the nucleus. Binding 
will occur in the special lowest angular momentum state J = / — | if 



''"■'^Tf 



(25) 



Thus binding can occur here only if the anomalous magnetic moment k > 1/4/. The proton, 
with K = 1.79, will bind. 

Binding can occur in higher angular momentum states J if and only if 



\k\ > Kc 



J^ + J-f 



(26) 



For example, for J = / + |, Kc = 2 + 3/4/, and for J = / + |, Kc = 4 + 15/4/. Thus ^He, 
which is spin 1/2, will bind in the first excited angular momentum state because k = —4.2. 
Unfortunately, to calculate the binding energy, one must regulate the potential at r = 0. 
The results shown in Table 1 assume a hard core. 



B. Nonrelativistic binding for general S 



The reference here is |^. The assumption made here is that / > S. (There are only 3 
exceptions, apparently: ^H, ®Li, and ^°B.) 

Binding in the lowest angular momentum state J = I — S is given by the same criterion 
(p5|) as in spin 1/2. Binding in the next state, with J = I — S + 1 occurs if A± > ;| where 



'-i)y-'^-'^^ 



[i + iy + i2s-i-i)^i + \pU 



(27) 



The previous result for S = 1/2 is recovered, of course. 5* = 1 is a special case: Then A_ is 
always negative, while A+ > ;j if 7 > 7c, where 



7 



^^(3 + 16/ + 16/^) 
^"4/9 + 4/ ^ ^ 



For higher spins, both X± can exceed 1/4: 



A+ > - for 7 > 7c„ (29) 

A_ > i for 7 > 7,+ (30) 



where for S* = | 



(7c)t = 4^(6 + 4/ T V33 + 32/). (31) 



For 4Be, for which 7 = —2.66, we cannot have binding because 3 > 7c_ > 1.557, 3 < 7c+ < 

5 

2' 



8.943, where the ranges come from considering different values of A^ from 1 to 00. For S — ^ 



, , 36 + 28/ + VII6I + 1296/ + 64/2 

(7c)^ = ^i • (32) 

So ^gAl will bind in either of these states, or the lowest angular momentum state, because 
7 = 7.56, and 1.67 > 7^- > 1.374, 1.67 < 7^+ < 4.216. 



C. Relativistic spin-1/2 



Kazama and Yang treated the Dirac equation p8|. See also ^9] and p3 



In addition to the bound states found nonrelativistically, deeply bound states, with 
-/^binding = Af are found. These states always exist for J > / + 1/2. For J = / — 1/2, 
these (relativistic) E = ^ bound states exist only if k > 0. Thus (modulo the question 
of form factors) Kazama and Yang ||28| expect that electrons can bind to monopoles. (We 



suspect that one must take the existence of these deeply bound states with a fair degree of 
skepticism. See also pO| .) 

As expected, for J = / — 1/2 we have weakly bound states only for k > 1/4/, which is 
the same as the nonrelativistic condition (pSf), and for J > / + 1/2, only if |fi;| > k^, where 
He is given in Eq. (^). 

D. Relativistic spin-1 

Olsen, Osland, and Wu considered this situation pT| . 

In this case, no bound states exist, unless an additional interaction is introduced (this is 
similar to what happens nonrelativistically, because of the bad behavior of the Hamiltonian 
at the origin). Bound states are found if an "induced magnetization" interaction (quadratic 
in the magnetic field) is introduced. Binding is then found for the lowest angular momentum 
state J = / — 1 again if k > 1/4/. For the higher angular momentum states, the situation is 
more complicated: 

8 



• for J = /: bound states require / > 16, and 

• for J > / + 1: bound states require J{J + 1) — /^ > 25. 

But these results are probably highly dependent on the form of the additional interaction. 
The binding energies found are inversely proportional to the strength A of this extra inter- 
action. 



Nucleus 


Spin 


7 


7 


J 


i?6 


Notes 


Ref 


n 


1 

2 


-1.91 




1 

2 


350 keV 


NR,hc 


m 


}H 


1 
2 


2.79 


2.79 


/-i = 


15.1 keV 

320 keV 

50-1000 keV 

263 keV 


NR,hc 

NR,hc 

NR,FF 

R 


m 
m 

INI 

i29| 


?H 


1 


0.857 


1.71 


/ - 1 = (iV = 2) 


If keV 


R,IM 


m 


|He 


1 
2 


-2.13 


-3.20 


'■ ^ 2 - 2 


13.4 keV 


NR,hc 


m 


f3^Al 


5 
2 


3.63 


7.56 


/-|=4 


2.6 MeV 


NR,FF 


m 


flAl 


5 
2 


3.63 


7.56 


/-|=4 


560 keV 


NR,hc 


133| 


ircd 


1 
2 


-0.62 


-1.46 


/+ 1 - 49 
'2 2 


6.3 keV 


NR,hc 


\n 



TABLE I. Weakly bound states of nuclei to a magnetic monopole. The angular momentum 
quantum number J of the lowest bound state is indicated. In Notes, NR means nonrelativistic and 
R relativistic calculations; he indicates an additional hard core interaction is assumed, while FF 
signifies use of a form factor. IM=induced magnetization, the additional interaction employed for 
the relativistic spin-1 calculation. We use A^ = 1 except for the deuteron, where A^ = 2 is required 
for binding. 



E. Remarks on binding 

Clearly, this summary indicates that the theory of monopole binding to nuclear mag- 
netic dipole moments is rather primitive. The angular momentum criteria for binding is 
straightforward; but in general (except for relativistic spin 1/2) additional interactions have 
to be inserted by hand to regulate the potential at r = 0. The results for binding energies 
clearly are very sensitive to the nature of that additional interaction. It cannot even be 
certain that binding occurs in the allowed states. In fact, however, it seems nearly certain 
that monopoles will bind to all nuclei, even, for example. Be, because the magnetic field in 
the vicinity of the monopole is so strong that the monopole will disrupt the nucleus and will 
bind to the nuclear, or even the subnuclear, constituents. 

F. Binding of monopole-nucleus complex to material lattice 

Now the question arises: Is the bound complex of nucleus and monopole rigidly attached 
to the crystalline lattice of the material? This is a simple tunneling situation. The decay 
rate is estimated by the WKB formula 



9 



r-V^rv/^M^^^, (33) 

a 

where the potential is crudely 

A4 is the nuclear mass <^ nionopole mass, and the inner and outer turning points, a and h 
are the zeroes oi E — V . Provided the following equality holds, 

{-Ef » ^, (35) 

which should be very well satisfied, since the right hand side equals 10~^^N^ MeV^, we can 
write the decay rate as 



r~ A^-i/^lO^V^exp 



3-137Vmp/ NB \m, 



(36) 



where the characteristic field, defined by cBq = m^, is 4 x 10 T. If we put in S = 1.5T, 
and A = 27, —E = 2.6MeV, appropriate for 13AI, we have for the exponent, for A^ = 1, 
—2 X 10^^, corresponding to a rather long time! To get a 10 yr lifetime, the binding energy 
would have to be only of the order of leV. Monopoles bound with kilovolt or more energies 
will stay around forever. 

Then the issue is whether the entire Al atom can be extracted with the 1.5 T magnetic 
field present in CDF. The answer seems to be unequivocally NO. The point is that the 
atoms are rigidly bound in a lattice, with no nearby site into which they can jump. A major 
disruption of the lattice would be required to dislodge the atoms, which would probably 
require kilovolts of energy [^. Some such disruption was made by the monopole when it 
came to rest and was bound in the material, but that disruption would be very unlikely to 
be in the direction of the accelerating magnetic field. Again, a simple Boltzmann argument 
shows that any effective binding slightly bigger than 1 eV will result in monopole trapping 
"forever." This argument applies equally well to binding of monopoles in ferromagnets. If 
monopoles bind strongly to nuclei there, they will not be extracted by 5 T fields, contrary to 
the arguments of Goto et al. |^ The corresponding limits on monopoles from ferromagnetic 
samples of Carrigan et al. |^ are suspect. 



IV. DUALITY AND THE EULER-HEISENBERG LAGRANGIAN 



Finally, let us consider the process contemplated in Refs. [0 and [|1J], that is 




+ 77, 77 -^ 77, (37) 

where the photon scattering process is given by the one-loop light-by-light scattering graph 
shown in Fig. |I|. If the particle in the loop is an ordinary electrically charged electron, this 

10 




FIG. 1. The light-by-light scattering graph for either an electron or a nionopole loop. 



process is well-known ||36|JT^j3^ . If, further, the photons involved are of very low momentuni 
compared the the mass of the electron, then the result may be simply derived from the well- 
known Euler-Heisenberg Lagrangian |3^, which for a spin 1/2 charged-particle loop in the 



presence of homogeneous electric and magnetic fields isQ 

1 /""^ ds __2 



C 



-T 



^TT^ Jo S^ 



• n^RecoshesX 2, ,,_ 

es) Q 1 (ts) T 

Im cosh esX 3 



(38) 



Here the invariant field strength combinations are 



T =\f- = \{^^ -^-) 



g= -F*F = E-H, 



(39) 



*F, 



fj-u 2 A"^Q/3 



F"^ being the dual field strength tensor, and the argument of the hyperbolic 
cosine in Eq. ( |38D is given in terms of 

X = [2(^ + iQ)]^/' = [(H + zE) Y^'. (40) 

If we pick out those terms quadratic, quartic and sextic in the field strengths, we obtain^ 

„2 1 



£ = -V+" 



4" ' 360 m^ 
vra^ 1 



[4(^2)2 ^ 'j^p*py 



630 m^ 



F^[8{Fy + 13(F*F)^] + 



(41) 



^We emphasize that Eq. (|3q ) is only valid when daF^y = 0. 
Incidentally, note that the coefficient of the last term is 36 times larger than that given in Ref. |1C| 
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The Lagrangian for a spin-0 and spin-1 charged particle in the loop is given by similar 
formulas which are derived in Ref. ^,16| and (implicitly) in Ref. |3^, respectively. 



Given this homogeneous-field effective Lagrangian, it is a simple matter to derive the 
cross section for the 77 — *> 77 process in the low energy limit. (These results can, of course, 
be directly calculated from the corresponding one-loop Feynman graph with on-mass-shell 



photons. See Refs. |ll6| , |37| .) Explicit results for the differential cross section are given by 
Ref. [§71: 



da 139 /I cu^ , 9 ^,n , , ^ 

-^ = ^«^ 3 + cos^ e?, (42) 

dn 324007r2 m^^ ' ' ^ ' 

and the total cross section for a spin-1/2 charged particle in the loop isQ 

973 ^u^ .._. 

a = a — -. (43 

101257r mS ^ ^ 

Here, uj is the energy of the photon in the center of mass frame, s = Auj"^. This result is 
valid provided u/m ^ 1. The dependence on m and u is evident from the Lagrangian (p|), 
the u dependence coming from the field strength tensor. Further note that perturbative 
quantum corrections are small, because they are of relative order 3a ~ 10~^ [0. Processes 
in which four final-state photons are produced, which may be easily calculated from the 
last displayed term in Eq. (^T]), are even smaller, being of relative order ~ a'^{uj/m)^. So 
light-by-light scattering, which has been indirectly observed through its contribution to the 
anomalous magnetic moment of the electron ^1|, is completely under control for electron 
loops. 

How is this applicable to photon scattering through a monopole loop? At first blush 
this calculation seems formidable. The interaction of a magnetically charged particle with 
a photon involves a "string," as described by the function /^ given in Eq. (|^). The inter- 
action between electric and magnetic charges is given by the complicated expression (^. 
This coupling is equivalent to the interaction between the magnetic current *J^ and the 
electromagnetic field, 

w^int = J{dx){dx'yF^,{x')r{x' - xyj^ix). (u) 

From Eqs. (|) and (^ one obtains the relevant string-dependent monopole-photon coupling 
vertex in momentum space, 

r,(g) = ^ /^-^'^^';^\ (45) 

n- q — le 

where we have, for variety's sake, chosen a semi-infinite string. As we have noted, the choice 
of the string is arbitrary; reorienting the string is a kind of gauge transformation. In fact, 
it is this requirement that leads to the quantization conditions (P and @). 



^The numerical coefficient in tlie total cross section for a spin-0 and spin-1 charged particle in the 
loop is 119/202507r and 2751/2507r, respectively. Numerically the coefficients are 0.00187, 0.0306, 
and 3.50 for spin 0, spin 1/2, and spin 1, respectively. 
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The authors of Refs. [0, |]T2[, and |T3| do not attempt a calculation of the "box" diagram 
with the interaction (|4^) . Rather, they (explicitly or implicitly) appeal to duality, that is, 
the symmetry that the introduction of magnetic charge brings to Maxwell's equations: 

E ^ H, H ^ -E, (46) 

and similarly for charges and currents. Thus the argument is that for low energy photon 
processes it suffices to compute the fermion loop graph in the presence of zero-energy pho- 
tons, that is, in the presence of static, constant fields. The box diagram shown in Fig. [^ with 
a spin-1/2 monopole running around the loop in the presence of a homogeneous E, H field 
is then obtained from that analogous process with an electron in the loop in the presence 
of a homogeneous H, — E field, with the substitution e —>■ g. Since the Euler-Heisenberg 
Lagrangian (^Tj) is invariant under the substitution (^61) on the fields alone, this means we 
obtain the low energy cross section cr^^^^^ through the monopole loop from Eq. ( P5| ) by the 
substitution e ^ g, or 

137 
a^ag = —-N\ N =1,2,3,.... (47) 



A. Inconsistency of the Duality Approximation 

It is critical to emphasize that the Euler-Heisenberg Lagrangian is an effective Lagrangian 
for calculations at the one fermion loop level for low energy, i.e., uj/M <^ 1. It is commonly 
asserted that the Euler-Heisenberg Lagrangian is an effective Lagrangian in the sense used 
in chiral perturbation theory [^,^]. This is not true. The QED expansion generates 



derivative terms which do not arise in the effective Lagrangian expansion of the Euler- 
Heisenberg Lagrangian [^. One can only say that the Euler-Heisenberg Lagrangian is a 
good approximation for light-by-light scattering (without monopoles) at low energy because 
radiative corrections are down by factors of a. However, it becomes unreliable if radiative 
corrections are large. 

In this regard, both the Ginzburg ||T2| , p!3[| and the De Riijula [l^ articles, particularly 



Ref. [0, are rather misleading as to the validity of the approximation sketched in the 
previous section. They state that the expansion parameter is not g but guj/M , M being the 
monopole mass, so that the perturbation expansion may be valid for large gf if a; is small 
enough. But this is an invalid argument. It is only when external photon lines are attached 
that extra factors of uj/M occur, due to the appearance of the field strength tensor in the 
Euler-Heisenberg Lagrangian. Moreover, the powers of g and uj/M are the same only for the 
F'^ process. The expansion parameter is Ug, which is huge. Instead of radiative corrections 
being of the order of a for the electron- loop process, these corrections will be of order Ug, 
which implies an uncontrollable sequence of corrections. For example, the internal radiative 
correction to the box diagram in Fig. ^ have been computed by Ritus |^ and by Renter, 



^The same has been noted in another context by Bordag et al. 
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Schmidt, and Schubert |^6| in QED. In the O(a^) term in Eq. (^) the coefficients of the 



(F2)2 and the (FF)^ terms are multiphed by (l + f f + ©(a^)^ ^nd (l + ^f + 0( 
respectively. The corrections become meaningless when we replace a ^ Ug. 



a^),, 



This would seem to be a devastating objection to the results quoted in Ref. |T3[ and used 
in Ref. fl^. But even if one closes one's eyes to higher order effects, it seems clear that the 
mass limits quoted are inconsistent. 

If we take the cross section given by Eq. (^31) and make the substitution (|47|) , we obtain 
for the low energy light-by- light scattering cross section in the presence of a monopole loop 



973 N^ uj^ . „ _, _« 1 



a. 



77^77 



25920007r a^ m^ M^ \M 



uo 



4.2xl0^iV«— ^ . (48) 



If the cross section were dominated by a single partial wave of angular momentum J, the 
cross section would be bounded by 

^7r(2J+l) Stt ,^^. 

a < — ~ — , (49) 

s s 

if we take J = 1 as a typical partial wave. Comparing this with the cross section given 



in Eq. (^8]) , we obtain the following inequality for the cross section to be consistent with 
unitarity, 

— > 3iV. (50) 



But the limits quoted [|T4| for the monopole mass are less than this: 

^ > 870 GeV, spin 1/2, (51) 

because, at best, a minimum {uj) ~ 300 GeV; the theory cannot sensibly be applied below 
a monopole mass of about 1 TeV. (Note that changing the value of J in the unitarity limits 
has very little effect on the bound (^Up since an 8th root is taken: replacing J by 50 reduces 
the limit (0) only by 50%.) 

Similar remarks can be directed toward the De Riijula limits [10|. That author, however. 



notes the "perilous use of a perturbative expansion in gP However, although he writes down 
the correct vertex, Eq. (^51) , he does not, in fact, use it, instead appealing to duality, and 
even so he admittedly omits enormous radiative corrections of 0{ag) without any justifica- 
tion other than what we believe is a specious reference to the use of effective Lagrangian 
techniques for these processes. 



B. Proposed Remedies 

Apparently, then, the formal small uj result obtained from the Euler-Heisenberg La- 
grangian cannot be valid beyond a photon energy uj/M ^0.1. The reader might ask why 
one cannot use duality to convert the monopole coupling with an arbitrary photon to the 
ordinary vector coupling. The answer is that little is thereby gained, because the coupling 
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of the photon to ordinary charged particles is then converted into a comphcated form anal- 
ogous to Eq. (44). This point is stated and then ignored in Ref. [1^] in the calculation of 
Z —^ 37. There is, in general, no way of avoiding the complication of including the string. 

We are currently undertaking realistic calculations of virtual (monopole loop) and real 
(monopole production) magnetic monopole processes. These calculations are, as the reader 
may infer, somewhat difficult and involve subtle issues of principle involving the string, and 
it will be some time before we have results to present. Therefore, here we wish to offer 
plausible qualitative considerations, which we believe suggest bounds that call into question 
the results of Ginzburg et al. |T^JTB| . 

Our point is very simple. The interaction ( PI) couples the magnetic current to the dual 
field strength. This corresponds to the velocity suppression in the interaction of magnetic 
fields with electrically charged particles, or to the velocity suppression in the interaction 
of electric fields with magnetically charged particles, as most simply seen in the magnetic 
analog of the Lorentz force. 



F = (7(B--xE). 



(52) 



That is, the force between an electric charge e and magnetic charge g^ moving with relative 
velocity v and with relative separation r is 



6(7 V X r 
c 47rr^ 



(53) 



This velocity suppression is refiected in nonrelativistic calculations. For example, the energy 
loss in matter of a magnetically charge particle is approximately obtained from that of a 
particle with charge Ze by the substitution [E^] 



V 



9_ 
c 



(54) 



And the classical nonrelativistic dyon-dyon scattering cross section near the forward direction 



is 19 



da 

'dn 



i2fivy 



ei92 - 62^1 

Attc 



+ 



r 6162 + 9192 

\ Atxv 



{e/2y 



<i, 



(55) 



the expected generalization of the Rutherford scattering cross section at small angles. 

Of course, the true structure of the magnetic interaction and the resulting scattering cross 
section is much more complicated. For example, classical electron-monopole or dyon-dyon 
scattering exhibits rainbows and glories, and the quantum scattering exhibits a complicated 
oscillatory behavior in the backward direction [^. These refiect the complexities of the 
magnetic interaction between electrically and magnetically charged particles, which can 
be represented as a kind of angular momentum [pO| , ^ . Nevertheless, for the purpose of 
extracting qualitative information, the naive substitution. 



-9^ 



(56) 
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seems a reasonable first step.[] Indeed, such a substitution was used in the proposal ^ to 
estimate production rates of monopoles at Fermilab. 

The situation is somewhat less clear for the virtual processes considered here. Never- 
theless, the interaction (0) suggests there should, in general, be a softening of the vertex. 
In the current absence of a valid calculational scheme, we will merely suggest two plausible 



alternatives to the mere replacement procedure adopted in Refs. [|T0],|12],|13|,|T5|. 

We first suggest, as seemingly Ref. |]T3[ does, that the approximate effective vertex in- 
corporates an additional factor of uj/M. Thus we propose the following estimate for the 77 
cross section in place of Eq. (ff3) , 



a. 



77^77 



rsj 



1 f UJ\ 
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'''^W[ff) • («^' 



since there are four suppression factors in the amplitude. Now a considerably larger value 
of u is consistent with unitarity, 

— > v^, (59) 

CO 

if we take J = 1 again. We now must re-examine the app^^^x cross section. 

In the model given in Ref. |TB|, where the photon energy distribution is given in terms 
of the functions f{y), y = uj/E, the physical cross section is given by 

(^pp-.-r-rX = {-) / V (yi)/(l/2)tT^7^T,T, = / dyidy2-: — 1— , (60) 

KtcJ J yi y2 J dyidy2 



where now (c/. Eq. (25) of in the first reference in Ref. ||T3[| ) 

da fa\'^^^f.fE\ 



dyidy2 Vvr 
where, for spin 1/2, (up to factors of order unity) 



) ^^' (m) 2/1^(^1)^2/(2/2), (61) 






This, and the extension of this idea to virtual processes, leaves aside the troublesome issue of 
radiative corrections. The hope is that an effective Lagrangian can be found by approximately 
integrating over the fermions which incorporates these effects. A first estimate of the effect of 
incorporating radiative correction, however, may be made by applying Pade summation of the 
leading corrections found in Refs. |45|,46|: 



^PT ^PT 



''^ {i-agY^ woom' ^^^^ 



which reduces the quoted mass Umits of Ref. ||lj] a bit. 
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The result in (|6TD differs from that in Ref. [|T3| by a factor of {E/M)^yfy2. The photon 
distribution function y'^f{y) used is rather strongly peaked at y ~ 0.3. (This peaking is 
necessary to have any chance of satisfying the low-frequency criterion.) When we multiply 
by 1/^, the amplitude is greatly reduced and the peak is shifted above y = 1/2, violating 
even the naive criterion for the validity of perturbation theory. Nevertheless, the integral of 
the distribution function is reduced by two orders of magnitude, that is, 

lodyy'fiy) ^ ' 



This reduces the mass limit quoted in [|T^] by a factor of 1/a/3, to about 500 GeV, where 
{uj)/M ^ 0.9. This dubious result makes us conclude that it is impossible to derive any 
limit for the monopole mass from the present data. 

As for the De Riijula limit^from the Z —>■ 3'-/ process, if we insert a suppression factor of 
co/M at each vertex and integrate over the final state photon distributions, given by Eq. (18) 
of Ref. |iy], the mass limit is reduced to M/^/N ^ \Amz ~ 120 GeV, again grossly violating 



the low energy criterion. And the limit deduced from the vacuum polarization correction to 
the anomalous magnetic moment of the muon due to virtual monopole pairs |I5| is reduced 
to about 2 GeV. 

The reader might object that this uj/M softening of the vertex has little field-theoretic 
basis. Therefore, we propose a second possibility that does have such a basis. The vertex 
(^) suggests, and detailed calculation supports (based on the tensor structure of the photon 

amplitudes^) the introduction of the string-dependent factor Jq^ /{n ■ qY at each vertex, 
where q is the photon momentum. Such a factor is devastating to the indirect monopole 
searches — for any process involving a real photon, such as that of the DO experiment |l^ 



or for Z ^ 37 discussed in |T^, the amplitude vanishes. Because such factors can and do 
appear in full monopole calculations, it is clearly premature to claim any limits based on 
virtual processes involving real final-state photons. 



^We note that De Riijula also considers the monopole vacuum polarization correction to gv IdA-, 
QA, and mw/iT^z, proportional to (mz/M)'^ in each case, once again ignoring both the string and 
the radiative correction problem. He assumes that the monopole is a heavy vector-like fermion, 
and obtains a limit of M/N > 8mz- Our ansatz changes {mz/M)"^ to {mz/M)'^, so that M/\/N > 
■\/8mz ~ 250 GeV, a substantial reduction. 

^For example, the naive monopole loop contribution to vacuum polarization differs from that of 
an electron loop (apart from charge and mass replacements) entirely by the replacement in the 
latter of (g^,y — q^q^/q'^) — > (q^/?o)('^ij ~ qiqj/^)^ when n^ points in the time direction. Apart 
from this different tensor structure, the vacuum polarization is given by exactly the usual formula, 
found, for example in Ref. |lq]. Details of this and related calculations will be given elsewhere. 
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V. CONCLUSIONS 

The field theory of magnetic charge is still in a rather primitive state. Indeed, it has 
been in an arrested state of development for the past two decades. With serious limits now 
being announced based on laboratory measurements, it is crucial that the theory be raised 
to a useful level. 

At the present time, we believe that theoretical estimates for the production of real 
monopoles are more reliable than are those for virtual processes. This is because, in effect, 
the former are dominated by tree-level processes. We have indicated why the indirect limits 
cannot be taken seriously at present; and of course only the real production processes offer 
the potential of discovery. Perhaps the arguments here will stimulate readers to contribute 
to the further development of the theory, for it remains an embarrassment that there is no 
well-defined quantum field theory of magnetic charge. 
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